Modeling electromagnetic form factors of light and heavy pseudoscalar mesons 
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The electromagnetic form factors of light and heavy pseudoscalar mesons are calculated within 
two covariant constituent-quark models, viz., a light-front and a dispersion relation approach. We 
investigate the details and physical origins of the model dependence of various hadronic observables: 
the weak decay constant, the charge radius and the elastic electromagnetic form factor. 



I. INTRODUCTION 

The light-front quantum field theory is discussed by 
N. Bogoliubov et al. in Ref. fH and an early applica- 
tion to study hadronic quark bound states may be found 
in Refs. [1, [3|. A self-consistent relativistic treatment 
of the quark spins can be performed within the light- 
front quark model, which allows for a calculation of the 
reference-frame independent partonic contribution to the 
form factor. The non-partonic contribution cannot be 
generally calculated, however it can be eliminated for 
space-like momentum transfers by an appropriate choice 
of reference frame. Therefore, the partonic contribution 
obtained in this specific reference frame yields the full 
form factor. 

The light-meson sector, which includes the pion and 
kaon, allows us to test QCD hypotheses on the subatomic 
structure of hadrons at low and intermediate energies. 
Yet, many aspects of quantum field theory on the light- 
front and its application to bound-state systems give rise 
to various open questions; for example, the problems 
of regularization and renormalization on the light-front. 
For the pion and kaon, the light-front constituent quark 
model (LFCQM) has been studied 0, i] to describe re- 
cent experimental data d, 0, H, Q • 

The description of these bound states using dispersion 
relations (DR) was first developed in JO] . It was applied 
to calculate light mesons form factors ll|, and used to 
elucidate long-distance effects [l^ [l4 1 in weak decays of 
heavy mesons. More recently the interplay between per- 
turbative and nonperturbative regions in the pion elec- 
tromagnetic form factor was discussed within the DR ap- 
proach [l^, where it was pointed out the dominance of 
the nonperturbative contribution up to rather high values 
of momentum transfers, going beyond a previous analy- 
sis performed up to intermediate values of momentum 
transfers within light-cone sum rule approach [iGj . 

The DR approach is based on a consistent treatment of 
the two-particle singularities which arise in the triangle 
diagrams describing elastic as well as inelastic meson- 
transition amplitudes. In this dispersive approach, these 
amplitudes are given by relativistic spectral integrals over 
the mass variables in terms of Bethe-Salpeter amplitudes 



of the mesons and spectral densities of the correspond- 
ing triangle diagrams. Thus, elastic form factors are de- 
scribed by double spectral representations. 

In this paper, we analyze the model dependence which, 
due to vertex functions as well as the choice of a con- 
stituent quark mass, usually arises in quark model cal- 
culations of elastic form factors. In the light sector, we 
compute the pion and kaon electric form factors for which 
experimental data is available; we then use our models 
to make predictions for the charmed sector. For sim- 
plicity we will assume point-like constituent quarks, as 
our model assumptions are made directly on the explicit 
form of wave functions and vertices. However, the con- 
stituent quark may have a structure that depends on the 
dynamical model of the bound state (see e.g. [l3l)- 



II. ELECTROMAGNETIC FORM FACTORS 

The space-like electromagnetic form factor of a pseu- 
doscalar meson with mass M is generally given by the 
covariant expression 



(1) 



with = p'^ = and the four-momentum transfer 
q = p' — p, < 0. The electromagnetic current is 
— 9(0)7/^9(0) where q{0) denotes a current quark. 
Here, F°™((7^) describes the virtual photon emission (ab- 
sorption) amplitude by the composite state of a quark 
of charge ei and an antiquark of charge 62- This form 
factor depends on the three independent Lorentz invari- 
ants p^, p'"^ and q^. The constituent quark amplitude of 
the electromagnetic form factor is assumed to have the 
following structure 

{Q{p')\miM0)\Qip)) = Qip')iMp)Uq') , (2) 

where Q{p) represents the constituent quark. The form 
factor is normalized such that F°™(0) = ei + 62 and 
we neglect the anomalous magnetic moment of the con- 
stituent quark in Eq. The function £,c{q^) describes a 
constituent quark transition form factor. Since the quark 
model is not formally derived from QCD, it is unknown. 
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In the following we assume that — Cc = ei or 62 owing 
to the fact that constituent quarks behave like bare Dirac 
particles [l8l |. 



III. LIGHT-FRONT CONSTITUENT QUARK 
MODEL 

The electromagnetic form factor for pseudoscalar par- 
ticles is calculated in the impulse approximation, using 
the Mandelstam formula 



eiiVc 



Tr[755i(fc-p') X 



(p + p')t^ J (27r)4 
rSi{k - ph^S^ik)] A{k,p)A{k,p') + [1^2], (3) 

where Si (p) is the Feynman propagator of quark i with 
constituent mass m^, Nc the number of colors, and 
A{k,p) is our hadron-quark vertex function model. The 
bracket [1 ^ 2] is a shortcut for the subprocess of con- 
stituent quark 2 interacting with the photon. 

The momentum component Bethe-Salpeter (BS) ver- 
tex model is chosen such that it regularizes the amplitude 
of the photo-absorption process and constructs a light- 
front valence wave function 0, 0, [l^. In the present 
study, we use two different models for the vertex func- 
tions taken from previous work, one of which has a non- 
symmetrical form under the exchange of the momentum 
of the quark and antiquark [3| , 



A(fc,p) = 



TV 



(p — k)^ — -|- ie ' 
while the other has a symmetric form 



A(fc,p) 



N 



N 



{p^ky 



- le 



(4) 



(5) 



The normalization constant N is obtained from (0) = 
ei -I- 62. 

The light-front constituent quark model for hadrons 
used here is based on quantum field theory and our ansatz 
comes with the choice of the BS vertex and point-like 
quarks. It is worthwhile to mention that the present 
LFCQM reproduces in its full complexity a covariant cal- 
culation. Technically, the use of the Drell-Yan frame (lon- 
gitudinal momentum transfer ~ q'^ + = 0) and the 
projection on to the light-front (performing the integra- 
tion on the minus momentum component analytically in 
the Mandelstam formula) simplifies drastically the com- 
putation of the form factor. This ap pro ach has been 
applied to the pion [1, [^, [l^ , the kaon ]20j as well as to 
the p 

The use of light-front variables in the evaluation of the 
Mandelstam formula, with the corresponding projection 
on the light-front hypersurface through the k^ — k'^ ~ k^ 
integration, has its subtleties. Covariance of the start- 
ing expression in Eq. ([3]), which corresponds to a frame 
independent form factor, is an important property to be 



checked in the final results. In some cases, as for vector 
mesons in the present model, it is necessary to perform a 
careful analysis of the light-front calculation of the photo- 
absorption amplitude which accounts for pair terms or Z- 
diagrams that survive in the Drell-Yan frame (see, e.g., 
Ref. [131). The Z-diagram vanishes in the calculation of 
the plus component of the electromagnetic current (j^) 
with — for pseudoscalar mesons with 76 coupling 
and the momentum component of the BS vertices from 
Eqs. Q and ©. 

Therefore, we adopt the plus-component convention of 
the electromagnetic current in the Breit-frame with the 
Drell-Yan condition (g"*" = 0) to compute the form factor 
of pseudoscalar mesons. The k~ integration is performed 
analytically by evaluating the residues in Eq. In 
the case of the non-symmetric vertex the resulting 
expression for the electromagnetic form factor is: 



i^™(g2) 



N^Nr 



ei 



dPk\dk+ Tr+[ 



{p+ +p'+) J 47r3 k+{p+ 



VP - k. 



2on {P ^)lon) {p' ^2on 



[1^2], 



{P' - k)ion) 

(6) 

[kl 



where the on-minus-shell momentum is k^^^ 
m\)/k'^ (^2on — (fcl + ^^i)/^^) ^''^d the Dirac trace from 
Eq. ^ for /i = + (good component of the current) is 



Tr^ 



-4 [{p+ - k+) (2mim2 + {p + p') - k- 2k^) 



-k+ {{p -k)-{p-k)- ml) 



(7) 



The electromagnetic form factor written for the non- 
symmetric vertex (|4]) in terms of the Bjorken momentum 
fraction x = k'^ jp^ is given by 



ei 



(Pk±dx 
87r3 J x{l - a;)4 

^{x,k±,qj_, nil, 1112) 



(M2 - M^){M'^ - M'2(m2, mR)) 
1 

[AP - M^){AP - M2(m2, mfl)) ^ 



[1^2], (8) 



where ^(x,fcj_,(f_L,"ii, TO2) = (p+) ^Tr^i ]\k-=k-^j ^ < 
x < 1 and the free two-quark mass is Mq = Af^(m2, ?ni). 



with 



if {p - k)\ + 



l-x 



^-pI 



(9) 



Analogous expressions follow for AI'^ and AI''^{m2,mii) 
with replaced by p'^. The calculation of the electro- 
magnetic form factor with the symmetrical vertex ([5]) can 
be done following the steps presented above. For the in- 
terested reader, the derivation of the form factor for the 
pion in this model is performed in detail in Ref. [5] . 

The mass parameters in these models are limited by 
TOi -|- 777.2 > AI, TOi -|- mR > M and TO2 + ''nR > AI 
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FIG. 1: The electromagnetic form factor in the impulse ap- 
proximation is obtained from the triangle diagram. The ver- 
tices and momenta depicted in the diagram represent the mo- 
mentum assignment in the DR approach. 



due to the unphysical continuum threshold as seen in the 
denominator of Eq. ([5]) . Quark confinement does not al- 
low for a scattering cut, but since in the present model 
the meson is as a real bound state, i.e. the conditions 
TOi -|- m2 > M and + m^j > M are satisfied, the cut is 
harmless. However, if the meson is a weakly bound com- 
posite particle, the form factor will be very sensitive to 
changes in the constituent quark mass. This is certainly 
not the case for pions and kaons, as they are Goldstone 
bosons and strongly bound, yet this sensitivity appears 
in the case of heavy pseudoscalars with small binding en- 
ergies in nonconfining models (see, e.g., Ref. [IS, H^])- 
particular, our numerical results for _D+ will exemplify 
the strong dependence of the heavy pseudoscalar electro- 
magnetic form factor on the constituent masses, whereas 



this effect is minor for confining models like the DR ap- 
proach. 



IV. DISPERSION RELATION APPROACH 

The amplitude is obtained from the triangle diagram 
shown in Fig. [1] where the kinematical variables are dis- 
played. The on-shell meson momenta are = and 
p'2 = M^. The triangle diagram may be calculated in 
various ways — we choose to put the constituent quarks 
on-mass shell while keeping the external momenta off- 
shell with 

p' = S, p'' = s', {p'-pf = q^. (10) 

Note, however, that p' — p = q ^ q. 

In the DR approach, the electromagnetic form factors 
^cmj-^2-j expressed by a double spectral representa- 
tion (see references [l^, Ell [13, El, E3|), 

' V 7r(s-M2) 7r(s'-M2) 

X A(s,s',92|mi,TO2) + [1^2]. (11) 

We apply the Landau-Cutkosky rules to calculate the 
double spectral density Ay: we place all internal parti- 
cles on their mass shell, k\ = m^, k'^ = mf, fc| = 
but take the variables and p'^ off-shell. The double 
spectral density is then derived from calculation of the 
triangle diagram 



2p^{q)A{s,s\q^\mi,m2) = ^ j d^kid^k[d''k2 8{kl - ml)5{kf 



ml)8{kl ~ ml) 



X 5(p - fci - k2)5{p' -k[- fca) Tr [-{jt[ + mi)7^(|^i + mi)i-f'^{m2 - h)^!^ 



(12) 



with (q) — Pfj, — ^ qf_i which ensures the Ward identity 
and thus charge conservation i^'^™(0) — ei + e2 with the 
proper vertex normalization. 

The vertex of the pseudoscalar meson in the con- 
stituent quark picture has the structure 



(5«(fci,mi)z75Q'^(-fc2,m2) 



G.(s) 



(13) 



where Q°'{ki,mi) represents the spinor state of the con- 
stituent quark of color a and iVj, the number of colors. 
The bound state vertex function may be related to the 
mesonic BS amplitude by 



(14) 



For a confining potential, the pole at s = Af^ should 



appear in the physical region for (j>{s) at s — Af^ > 
(mi -f 7712)^. This is, however, not the well known 

from the behavior of the bound-state wave function in an 
harmonic oscillator potential, 0(s) is a smooth exponen- 
tial function of s > (mi -I- m2)^. This means that the 
would-be pole in 0(s) at s = M'^ is completely blurred 
out by the interaction and it is therefore more appropri- 
ate to analyze the meson form factors in terms of (f){s) 
rather than G„(s). 

For a pseudoscalar meson in the dispersion approach 
of the constituent quark model, the BS amplitude cjj^s) 
which accounts for soft constituent quark rescattering is 
given by E3 




- (ml 



,2\2 



w{k) 
(mi — 7712)^ s'^/^ 



(15) 
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The dynamical factor multiplying 'w{k) stems from the 
loop diagram associated with the meson-vertex normal- 
ization. The modulus of the center of mass momentum 
is 



k = 



(s + mf — m^)^ — Asjnf 
4~s 



(16) 



It can be shown |13| that the vertex normalization of 
Gv{s) = 4>{s){s — M^), which describes soft constituent 
rescattering, reduces to a simple normalization of the 
wave function 



(17) 



A heuristic choice must be made for wlk). For phenom- 
ena that are predominantly governed by infrared mass 
scales, it is sensible in the case of heavy mesons to choose 
a BS amplitude parameterized by a function whose sup- 
port is in the infrared; for instance functions of Gaussian 
form, 



w{k) = Nexp [-Avk^/^?) 



(18) 



where the reduced mass of the quark-antiquark pair is 
/i = toiTO2/(toi -t- m2) and N a normalization constant. 
For light pseudoscalar mesons, dynamical chiral symme- 
try and the Ward-Takahashi identity are of great im- 
portance for their structure and properties. In non- 
perturbative approaches, using for example Schwinger- 
Dyson equations which incorporate these QCD features, 
appropriate parameterizations for the pion or kaon BS 
amplitudes are derived [1^. In this work, we employ 
power-law wave functions of the form 



Wn{k) 



N 



{l + v{k|^if 



n>2 



(19) 



inspired by the form suggested in Ref. [2y| from the anal- 
ysis of the ultraviolet physics of QCD. 

The size parameter u is to be determined from exper- 
imental and theoretical considerations. On the experi- 
mental side, on can either choose to constrain the vertex 
functions, or equivalently the BS amplitude, by the weak 
decay constant or by the electric charge radius. It is 
worthwhile to observe that both wave function models, 
Gaussian and power-law, do not present the singularity 
problem brought by the scattering cut due to changes the 
constituent quark masses, differently from the present 
LFCQM. 



V. NUMERICAL RESULTS 

The electromagnetic form factor of light and heavy 
pseudoscalar mesons are calculated with the two covari- 
ant constituent quark models. In the LFCQM, the quarks 
are in a bound state while in the DR model the wave 
functions corresponds to confined quarks with Gaussian 
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FIG. 2: The pion electormagnetic form factor in the LFCQM 
and DR approaches compared with experimental data from 
Ref. [^, Q, [3, 1^. The model and respective parameters are 
explained in the legend. 



and power-law forms. The shape of these functions is ob- 
tained from a fit to the experimental values of /^r , /k and 
fo or the charge radius (when available) of the regulator 
mass mR in the LFCQM and size parameter v in the DR 
approach. The light constituent quark masses uiu = rud 
are allowed to take the values 0.22 and 0.25 GeV (and 
0.28 GeV in case of the D^). The strange and charm 
constituent masses are fixed to = 0.508 GeV and 
rric = 1.623 GeV, respectively (see Ref. 

The results for the pion electromagnetic form factor 
compared to experimental data @, 0, @| are shown in 
Fig. [21 The quark masses used in these calculations are 
mu^md^ 0.22 GeV in both the LFCQM and DR ap- 
proaches. The pion charge radius is used to fit the param- 
eters TOfl and V (power-law wave function with n — Z). 
The regulator mass is m^j — 0.946 GeV for the non- 
symmetric vertex and 0.546 GeV 0] the symmetric ver- 
tex. In this case, we show results only for the power- 
law wave function with n — i {v = 0.088) in the DR 
approach because the Gaussian vertex produces a form 
factor strongly damped with increasing , which vio- 
lates the asymptotic QCD result predicting a q^'^ falloff. 
The pion electromagnetic form factor for the power-law 
wave function and non-symmetric vertex are very simi- 
lar, while the symmetrical model has a longer tail. These 
results indicate that the electromagnetic form factor of 
the pion as a strongly bound system does not particularly 
distinguish between models with a scattering cut. The 
high momentum experimental data appear to favor the 
more complex structure of the symmetrical vertex. It is 
expected that in this region the details of the regulator 
structure are crucial, and moreover a symmetric form of 
the vertex is natural. 

We remind that the tail of the momentum component 
of the pion light-front wave function corresponding to 
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FIG. 3: The kaon electromagnetic form factor in the LFCQM 
and DR a ppr oaches compared with experimental data from 
Refs. [2^. l29l|. The models and respective parameters are 
explained in the legend 

the non-symmetric 4] and to the symmetric vertex Q 
both decrease as ^ kj^, which can also be verified by 
inspecting Eq. ([8]). The power-law wave function model 
with n = 3 also displays a tail that falls off with ~ kj^. 
Thus, in the space- like region up to 10 GeV^, we find that 
the pion electromagnetic form factor favors the power-law 
damping of the wave function with a ^ fcj^ tail. 

The electromagnetic form factor of the obtained 
in both models are compared with experimental data 
psL I29I I in Fig. [3l The quark masses used in these calcu- 
lations are m„ — irid = 0.22 GeV and iris — 0-508 GeV 
in both the non-symmetric vertex and DR approaches. 
In order to study the mass dependence of the form fac- 
tor, we also give results for a calculation with m„ — 
TUd = 0.25 GeV in LFCQM case. The regulator mass 
is fixed to the pion's value mn = 0.946 GeV for the 
non-symmetric vertex which corresponds to a reasonable 
kaon charge radius [lOj. We observe in Fig. [3] that the 
change of 0.03 GeV in the constituent quark mass in is 
not important as it is much smaller than the kaon bind- 
ing energy of about 0.2 GeV. A similar finding can be 
reported for the DR approach, where the size parameter 
fi of the power-law wave function with n — 2 is fitted 
= 0.061) to reproduce the kaon charge radius and is 
readjusted for n = 3 (z^ = 0.158). Again, we observe that 
the non-symmetric vertex as well as the power-law model 
essentially yield the same form factor once the radius is 
reproduced. 

Our numerical results for Z?+ are shown in Fig. |4] for 
the DR approach and LFCQM. The size parameter of the 
Gaussian and power-law (n = 3) models are fitted to fo 
for irid = 0.25 GeV. Up to 10 GeV^ there is no sizable 
difference between the two calculations. This is reason- 
able as rric sets one large scale whereas the momentum 



FIG. 4: Predictions for the D electromagnetic form factor 
in the LFCQM and DR approach. 

transfer goes only up to 3 GeV, which is not large enough 
to discriminate between the models. The regulator mass 
of 1.77 GeV in the non-symmetric vertex model is cho- 
sen so as to give a sample of the LFCQM results. Here, 
we want to illustrate the constituent quark mass sensi- 
tivity when the symmetric and non-symmetric vertices 
are used. We remind that the mass parameters in these 
models are limited by mi -I- TO2 > M, mi + uir > M and 
"^2 + niR > M due to the unphysical continuum thresh- 
old as evident from the denominator in Eq. ([8]). The 
model binding energy of is only a few tenths of MeV, 
as can be inferred from the I?+ and quark masses. Since 
the meson mass is fixed to the experimental value of 
1.8694 GeV and rric = 1.623 GeV, the light quark mass 
rrid must be larger than the value 0.22 GeV we used for 
the pion and kaon. The increase of the light quark con- 
stituent mass from 0.25 to 0.28 GeV corresponds to a 
large increase in binding and consequently to a decrease 
of the radius, as seen in Fig. 31 

VI. CONCLUSION 

We have concentrated on the model dependence of elas- 
tic electromagnetic form factors of light and heavy pseu- 
doscalar mesons. Two different cases were studied by 
changing the form of the vertex or wave function and by 
varying the constituent quark masses. In the light sector, 
we computed the pion and kaon electric form factors for 
which experimental data is available; we then adjusted 
our models to make predictions for the charmed sector. 
Two ansatzes have been studied in the LFCQM which 
involve a symmetric and a non-symmetric vertex form. 
In the DR approach, Gaussian and power-law wave func- 
tions were used in the calculations. 

The form factors of the pion and kaon, both being 
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strongly bound systems, are sensitive to the short-range 
part of the vertex functions. Therefore, the tail of the 
momentum component of the light-front wave functions 
for the non-symmetric [3] and symmetric vertex 5] as 
well as the power-law model (n — 3), all decreasing as 
~ kj^, is essential to reproduce the pion space- like data 
up to 10 GeV^. Moreover, as a consequence of the strong 
binding of the quarks, changes in the constituent mass 
of about ten percent are not important for the elastic 
form factor of light mesons, provided the charge radius 
(or decay constant) is accordingly refitted. We have ex- 
emplified this in the kaon calculation. The light quark 
constituent masses are relevant for low momentum trans- 
fers of the order ^ ^ around 0.05 — O.SGeV^, which is 
dominated by the charge radius and consequently com- 
pensated by fitting the size or regulator mass parame- 
ters. Therefore, in qualitative agreement with QCD scal- 
ing laws in the ultraviolet limit, a power-law form of the 
wave function and our choice of vertices are reasonable in 
reproducing the form factor data once the charge radius 
is fitted. 

A recent calculation of the pion form factor within the 
dispersion relation approach including 0{as) QCD per- 
turbative contribution [15i] . showed in a model indepen- 
dent way that the nonperturbative part stays above 50% 
for < 20 GeV^ . The nonperturbative part is compati- 
ble with a pion wave function close to the asymptotic one, 
with a momentum scale given by the effective continuum 
threshold that determines to a great extent their results 
for the form factor. The value of the threshold is ex- 
tracted from the experimental pion decay constant. Simi- 
larly in our analysis the size or regulator mass parameters 
are fixed by the experimental values of the weak decay 
constants within models that have momentum tails de- 
creasing with ^ kj^, which does not leave not too much 
freedom for the pion form factor results below 20 GeV^, 
as we have shown. In that respect our calculation and 
[m show the dominance of the scale given by the pion 
decay constant that dials the vertex or wave function (de- 
caying as a power law) in the pion form factor up to a 
fairly large value of momentum transfers. 

The electromagnetic form factor of pseudoscalar 
mesons in the heavy-light sector is strongly sensitive to 
the infrared behavior of QGD, compared to the case of 
light mesons. The relevance of soft physics in this sec- 
tor is supported by our calculations of the elastic D me- 
son form factor both within the dispersion relation and 
light-front approaches. Once the size parameter of the 
Gaussian and power-law wave functions are adjusted to 
reproduce the decay constant, the elastic form factor is 



insensitive to the vertex model in the space-like region 
up to 10 GeV^. Therefore, the large- momentum tail of 
the wave function is not important for these results. We 
presume that the asymptotic behavior will be important 
only for ^ beyond the results we presented. The 
importance of soft physics may also be appreciated in the 
covariant calculations of the form factor with the sym- 
metric and non-symmetric vertices in LFCQM. The cor- 
responding BS amplitude models a lightly bound heavy- 
light system with respect to the light quark mass. The 
infrared dynamics appears to be important for the D 
meson, as the sensitivity to changes of the constituent 
quark mass by about 10% suggests. This small modifi- 
cation produces a sensible change in the nearby unphys- 
ical continuum threshold, as the mass of the D is about 
the sum of the constituent quark masses. Therefore, the 
condition mi -|- m2 > M is barely satisfied (see, e.g., the 
denominator of Eq. 

In summary, our study of the model dependence of 
elastic form factors for pseudoscalar mesons in the light 
and heavy-light sectors suggests a separation of the ul- 
traviolet and infrared physics. The ultraviolet properties 
of QGD dominate the form factor of light pseudoscalars, 
whereas infrared physics and details of quark confinement 
appear to be important for the space-like form factor of 
heavy mesons below = 10 GeV^. Hence, the electro- 
magnetic form factors of heavy-light systems, such as the 
D and B mesons, provide a valuable tool in the effort to 
investigate the nonperturbative physics of confinement, 
in contrast with light pseudoscalar form factors insensi- 
tive to the details of infrared physics even at low momen- 
tum transfers once the weak decay constants are fixed to 
the experimental values. 
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